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ABSTRACT 

We present a new technique for measuring the abundance of europium, a representative r-process 
element, in solar-metallicity stars. Our algorithm compares LTE synthetic spectra with high-resolution 
observational spectra using a ^-minimization routine. The analysis is fully automated, and therefore 
allows consistent measurement of blended lines even across very large stellar samples. We compare 
our results with literature europium abundance measurements and find them to be consistent; we also 
find our method generates smaller errors. 
Subject headings: data analysis and techniques — stars 



1. INTRODUCTION 

Every atom heavier than lithium has been processed 
by stars. Elements of the a process, iron peak, and r 
process have different formation sites, and therefore, un- 
derstanding the distribution of these elements in nearby 
stars can lead to a better understanding of the Galaxy's 
chemical enrichment history. The r-process site is the 
least well understood. Europium is our choice for an 
r-process investigation for two reasons: 96 % of Galactic 
europ ium is formed through the r-process (|Burris et al.l 
I2000T ). and it has several strong lines in the visible por- 
tion of the electromagnetic spectrum. In order to provide 
insight into the Galactic enrichment history, however, eu- 
ropium measurements in very large stellar samples will 
be needed. 

The study of r-process eleme nts is well estab- 
lished in metal -p oor stars (e . g.. |Sneden et alj 120081; 
Lai et all l2008t iFrebel et all 120071 : I Simmer er et all 
2004 Uohnson fc Boltd |2001[ ). where r-process abun- 
dances probe early Galactic history. A handful of 
studies have measured eu ropiu m in solar-meta llicity 
stars (e.g.. iReddv et all l2006t IBensbv et all 120051: 
iKoch fc Edvardssonll2002TlWoolf et al.lll995l) . but have 
done so in samples of 50-200 stars. We aim to extend 
such work to a substantial sample of 1000 stars at solar 
metallicity, using spectra collected by the California and 
Caltech Planet Search (CCPS). Such a large sample will 
require automated analysis. 

Our goal in this paper is to establish a new, automated 
abundance f itting method based on S pectroscopy Made 
Easy (SME; IValenti fc Piskunovlll996D . an LTE spectral 
synthesis code discussed in further detail in jj[3l The 
analy sis builds on the framework of the lValenti fc Fischerl 
120051 (hereafter VF05) Spectroscopic Properties of Cool 
Stars (SPOCS) catalog. Our method will yield consistent 
results across large stellar samples, generating smaller 
errors than previous analyses. 

We begin here with 41 stars that have been examined 
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in the literature, and we include three europium lines: 
4129 A, 4205 A, and 6645 A. Our stellar observations are 
detailed in §0 We describe the details of our europium 
abundance measurement algorithm in §[3l compare our 
values with existing europium literature in §|4j and sum- 
marize the results of the study in §[5] 

2. OBSERVATIONS 

The spectra of the 41 stars included in this study were 
taken with the HIRES echelle spectrograph (jVogt et al.l 
1994) on the 10-m Keck I telescope in Hawaii. The spec- 
tra date from January 2004 to September 2008 and have 
resolution R ~ 50 000 and signal-to-noise ratio (S/N) of 
~ 160 near 4200 A, where the two strongest Eu II lines 
included in this study are located. The spectra have the 
same resolution but S/N ~ 340 at 6645 A, where a third, 
weaker Eu n line is located. 

The spectra were originally obtained by the CCPS with 
the intention of detecting exoplanets. The CCPS uses 
the same spectrometer align ment each night and employs 
the HIRES exposure meter (jKibrick et al.ll2006h . so the 
stellar observations are extremely consistent, even across 
years of data collection. For a mo re complete descr iption 
of the CCPS and its goals, see iMarcy et all 120081 It 
should be noted that t he iodine cell Doppler technique 
(|Marcv fc Butlerl I1992D imprints molecular iodine lines 
only between the wavelengths of 5000 and 6400 A, leaving 
the regions of interest for this study iodine free. 

2.1. Stellar Sample 

For this study we choose to focus on a subset of 
stars that have been analyzed in other abundance stud- 
ies. We compare our measurements wi t h those from 
Woolf et al.l 1 19951 ISimmerer et al.l 120041 : IBensbv et all 
20051: and ldel Peloso et al.ll2005l We select CCPS target 
stars that have temperature, gravity, mass, and metal- 
licity measurements in the SPOCS catalog (VF05). 

Our initial stellar sample consisted of 44 objects that 
were members both of the aforementioned studies and 
the SPOCS catalog, and that were observed by the 
CCPS on the Keck I telescope. Three stars that oth- 
erwise fit our criteria, but for which we could only mea- 
sure the europium abundance in one line, were removed 
from the sample. They were HD 64090, HD 188510, and 
HIP 89215. While stellar europium abundances may be 
successfully measured from a single line, the goal of this 
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study is to establish the robustness of our method, and 
so we deem it necessary to determine the europium abun- 
dance in at least two lines to include a star in our analysis 
here. We describe our criteria for rejecting a fit in § 13.21 
The 41 stars included in this study have 3.4 < V < 
10.0, 0. 50 < B - V < 0.92, and 5 < d < 75 pc (ffip- 
parcos; lESAl [19971. VF05 determine stellar properties 
using SME, so it is reasonable to adopt their values for 
our SME europium analysis. Based on the VF05 anal- 
ysis, our stars have metallicity —1.4 < [M/H] < 0.4, 
effective temperature 4940 < T c g < 6230, and gravity 
3.8 < log g < 4.8. It should be noted that throughout 
this study we use the [M/H] parameter for a star's metal- 
licity, rather than the iron abundance [Fe/H]. Our [M/H] 
value is taken from VF05, where it is an independent 
model parameter that adjusts the relative abundance of 
all metals together. It is not an average of individual 
metal abundances. The full list of stellar properties, from 
the Hipparcos and SPOCS catalogs, appears in Table Q] 

2.2. Co-adding Data 

The nature of the radial velocity planet search dictates 
that most stars will have multiple, and in some cases, 
dozens of, observations. To take advantage of the mul- 
tiple exposures, we carefully co-add the reduced echcllc 
spectra where possible. 

The co- adding procedure is as follows: a 2000- pixel re- 
gion (approximately half the full order width) near the 
middle of each order is cross-correlated, order by order, 
with one observation arbitrarily designated as standard. 
The pixel shifts are then examined and a linear trend 
as a function of order number is fit to the pixel shifts. 
For any order whose pixel shift falls more than 0.4 pixels 
from the linear trend, the value predicted by the linear 
trend is substituted. This method corrects outlying pixel 
shift values, which often proved to be one of a handful 
of problematic orders where the echelle blaze function 
shape created a false cross-correlation peak. 

Each spectral order is adjusted by its appropriate frac- 
tional pixel shift before all the newly aligned spectra are 
added together. In order to accurately add spectra that 
have been shifted by non-integer pixel amounts, we use 
a linear interpolation between pixel values to artificially 
(and temporarily) increase the sampling density by a fac- 
tor of 20. After the co-adding, the resultant high-S/N 
multi-observation spectrum is sampled back down to its 
original spacing. 

The number of observations used per star is recorded in 
the N b s column of Table [2J The co-adding proved par- 
ticularly beneficial when fitting the relatively weak Eu II 
line at 6645 A. Two sample stellar spectra, from a star 
with 1 observation and from a star with 17 observations, 
are given in Figure [1] as evidence of the S/N advantages 
of this procedure. 

3. ABUNDANCE MEASUREMENTS 

We use the SME suite of routines for our spectral syn- 
thesis, both for fine-tuning line lists based on the solar 
spectrum (§ 13. lj) and for measuring europium in each star 
(§ 13.21). SME is an LTE spectral synthesis code based on 
the lKuruczlll992l grid of stellar atmosphere models. In 
brief, to produce a synthetic spectrum, SME interpolates 
between the atmosphere models, calculates the continu- 
ous opacity, computes the radiative transfer, and then 
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Fig. 1. — The three Eu II lines considered in this paper, plot- 
ted for two different stars to demonstrate the advantage to be 
gained from co-adding multiple observations of the same star. 
Note that the rightmost plots have a different ordinate axis scal- 
ing than the other four panels. The top three panels, showing 
HD 156365, include 1 spectrum, while the bottom three panels, 
showing HD 144585, include 17. The two stars have approximately 
the same metallicity and effective temperature. The advantage 
from co-adding is most profound in the weak 6645-A line. 

applies line broadening, which is governed by macrotur- 
bule nce, stellar rotation, and instrumental profile. Con- 
sult Vale nti fc Piskuno v 1996 and VF05 for a more m- 
depth description of SME's inner workings. 

Throughout this study we use SME only to compute 
synthetic spectra. All fitting is done in specialized rou- 
tines of our own design, external to SME. 

In this section we outline our technique for measur- 
ing europium abundances in the set of 41 stars included 
in this work. In broad strokes, we first determine the 
atomic parameters of our spectral lines by fitting the 
solar spectrum (§ 13. ip . Then, we use that line list to 
measure the europium abundance in three transitions 
(4129 A, 4205 A, and 6645 A); a weighted average of the 
three transitions determines a star's final europium value 
(§ 13. 2p . Finally, we estimate our uncertainties by adding 
artificial noise to our data in a series of Monte Carlo tri- 
als f§ !3.3p . We also here include notes on individual lines 

(MM- 

3.1. Line Lists 

We use relatively broad spectral segments in our eu- 
ropium analysis. The regions centered on the Eu II 
4129 A and 6645 A lines are 5 A wide, and the region 
centered on the Eu n 4205 A line is 8 A wide. We find it 
necessary to use such broad spectral segments in order 
to fit a robust and consistent continuum in the crowded 
blue regions. 

Line lists are initially drawn f rom the Vienna As- 
trophysics Line D atabase (VALD: IPiskunov et al.l fl995h 
iKupka et aLlll999D . The VALD line lists, in the regions 
surrounding all three Eu n transitions, make extensive 
use of Kurucz line lists. 

We apply the original VALD line list to an observed 
solar spectrum in order to determine the list's com- 
pleteness and to adjust line parameters as needed. In 
the blue, we use t he disk-center solar spectrum from 
IWallace et aLlll998l with the following global parameters: 
T eff = 5770 K, logg = 4.44, [M/H] = 0, microturbulence 
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TABLE 1 
Stellar Data. 



Identification 




B - V a 




T efI b 


f TV !f 1 T Tl Y~l 

[M/HJ 


log 9 b 


Ref. 


HD 


HR 


HIP 


(mag) 


(mag) 


(PC) 


(K) 


(cgs) 


3795 


173 


3185 


6.14 


0.72 


28.6 


5369 


-0.41 


4.16 


1 


4614 


219 


3821 


3.46 


0.59 


6.0 


5941 


-0.17 


4.44 


1 


6734 




5315 


6.44 


0.85 


46.4 


5067 


-0.28 


3.81 


1 


9562 


448 


7276 


5.75 


0.64 


29.7 


5939 


0.19 


4.13 


1, 2, 4 


9826 


458 


7513 


4.10 


0.54 


13.5 


6213 


0.12 


4.25 


4 


14412 


683 


10798 


6.33 


0.72 


12.7 


5374 


-0.45 


4.69 


1 


15335 


720 


11548 


5.89 


0.59 


30.8 


5891 


-0.20 


4.07 


A 


16397 




12306 


7.36 


0.58 


35.9 


5788 


-0.35 


4.50 


1 


22879 




17147 


6.68 


0.55 


24.3 


5688 


-0.76 


4.41 


1, 2, 4 


23249 


1136 


17378 


3.52 


0.92 


9.0 


5095 


0.03 


3.98 


1 


23439 




17666 


7.67 


0.80 


24.5 


5070 


-0.73 


4.71 


3 


30649 




22596 


6.94 


0.59 


29.9 


5778 


-0.33 


4.44 


1 


34411 


1729 


24813 


4.69 


0.63 


12.6 


5911 


0.09 


4.37 


4 


43947 




30067 


6.61 


0.56 


27.5 


5933 


-0.28 


4.37 


2 


45184 


2318 


30503 


6.37 


0.63 


22.0 


5810 


0.03 


4.37 


1 


48938 


2493 


32322 


6.43 


0.55 


26.6 


5937 


-0.39 


4.31 


1 


84737 


3881 


48113 


5.08 


0.62 


18.4 


5960 


0.14 


4.24 


1 


86728 


3951 


49081 


5.37 


0.68 


14.9 


5700 


0.11 


4.29 


4 


102365 


4523 


57443 


4.89 


0.66 


9.2 


5630 


-0.26 


4.57 


2 






57450 


9.91 


0.58 


73.5 


5272 


-1.42 


4.30 


3 


103095 


4550 


57939 


6.42 


0.75 


9.2 


4950 


-1.16 


4.65 


3 


109358 


4785 


61317 


4.24 


0.59 


8.4 


5930 


-0.10 


4.44 


4 


115617 


5019 


64924 


4.74 


0.71 


8.5 


5571 


0.09 


4.47 


4 


131117 


5542 


72772 


6.30 


0.61 


40.0 


5973 


0.10 


4.06 


2 


144585 




78955 


6.32 


0.66 


28.9 


5854 


0.25 


4.33 


1 


156365 




84636 


6.59 


0.65 


47.2 


5856 


0.24 


4.09 


1 


157214 


6458 


84862 


5.38 


0.62 


14.4 


5697 


-0.15 


4.50 


4 


157347 


6465 


85042 


6.28 


0.68 


19.5 


5714 


0.03 


4.50 


1 


166435 




88945 


6.84 


0.63 


25.2 


5843 


0.01 


4.44 


1 


169830 


6907 


90485 


5.90 


0.52 


36.3 


6221 


0.08 


4.06 


1 


172051 


6998 


91438 


5.85 


0.67 


13.0 


5564 


-0.24 


4.50 


1 


176377 




93185 


6.80 


0.61 


23.4 


5788 


-0.23 


4.40 


1 


179949 


7291 


94645 


6.25 


0.55 


27.0 


6168 


0.11 


4.34 


1 


182572 


7373 


95447 


5.17 


0.76 


15.1 


5656 


0.36 


4.32 


1 


190360 


7670 


98767 


5.73 


0.75 


15.9 


5552 


0.19 


4.38 


1 


193901 




100568 


8.65 


0.55 


43.7 


5408 


-1.19 


4.14 


3 


199960 


8041 


103682 


6.21 


0.64 


26.5 


5962 


0.24 


4.31 


1 


210277 




109378 


6.54 


0.77 


21.3 


5555 


0.20 


4.49 


1 


217107 


8734 


113421 


6.17 


0.74 


19.7 


5704 


0.27 


4.54 


1 


222368 


8969 


116771 


4.13 


0.51 


13.8 


6204 


-0.08 


4.18 


1 


224383 




118115 


7.89 


0.64 


47.7 


5754 


-0.06 


4.31 


1 



a ^-magnitude, color index, and parallax-based distance from the Hipparcos catalogue 
llESAI [1997), b Stellar parameters prev iously published in VF 05. C Star included for 
co mparison to the follow ing works: 1 — Bcnsbv ct al. 2005; 2 — del Pcloso ct all 120051 : 
3- JSimmerer et al.ll200l 4— IWoolf et al.1119951 



v-mic = 1.0 km s -1 , macroturbulence v mac = 3.6 km s _1 , 
rotational velocity vsmi = km s _1 , and radial velocity 
v r ad = km s _1 . These are the same sol ar parameters 
adop ted in VF05. In the red, we find the IWallace et ail 
119981 solar atlas to have insufficient S/N to accurately de- 
termine the atomic parameters. At 6645 A, therefore, we 
instead compare our line list to the disk-integrated NSO 
solar spectrum (|Kurucz et al.lll984l ). adjusting wsini to 
1.6 km s _1 because the full solar disk has more substan- 
tial rotational broadening. 

We find that adjustments to the oscillator strength 
(logg/) and van der Waals broadening (T 6 ) parameters 
are required for the strongest lines in a given wavelength 
segment, even far from the Eu II line of interest. For 
example, the Fe I line at 4202 A has an equivalent width 
W = 326 mA in the Sun and significantly affects the 
continuum of the 4205-A region. The log gf parameter 
controls the line depth while the Tg parameter controls 
the line shape, so in general the two parameters are or- 
thogonal. 



We used the Kurucz log gf values provided by VALD 
where possible, but adjustments were necessary where 
line depths were po orly fit. For Tg, VALD returns the 
iBarklem et al.l 120001 parameters for beryllium through 
barium (Z of 4-56), but has no T 6 values above Z — 56. 
We therefore find it necessary to fit Tq in species heavier 
than barium and in deep features not fit well by VALD 
values. We take particular care to determine the appro- 
priate log gf and Tg parameters for all lines adjacent to 
the Eu II line of interest. We find the best value for 
these parameters with the SME synthesizer by perform- 
ing a x 2 minimization against the solar spectrum on each 
parameter. 

Where the VALD line list is insufficient, we add line 
data from lMoore et al1 ll966. NIST (whose lists are based 
on a variety of sources) and C. Sneden (private com- 
munication). We also add CH and CN molecular lines 
based on values obtained from the Kurucz molecular line 
list web site. 3 In the 4129-A region, we find it neces- 

3 |http : //kurucz . harvard . edu/LIHELISTS/ LINESMOL/ 
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TABLE 2 
Stellar abundance values. 
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All names are HD numbers unless otherwise indicated. 15 For 



- loge(Eu) , 
of each line 



where 
in the 



a given star, [Eu/H] = loge(Eu) - 
logepO =log w (N x /N„). c The weight 
average is based on 50 Monte Carlo trials in each Eu II line. We 
have adopted an error floor of 0.02 dex, added in quadrat ure t o the 
errors determined by our Monte Carlo procedure. See i|3,3l for a 
more complete description of the weighted average and associated 
uncertainty. 



sary to add three artificial iron lines in order to match 
the solar spectrum. We follow here the precedent of 
Idel Peloso et al.ll2005l . though we find our fit requires the 
lines to have slightly different wavelengths and log gf val- 
ues. The complete line list for Eu 11 at 4129 A appears in 
Table H for 4205 A in Table gj and for 6645 A in Table 
03 The corresponding plots of the regions in these tables 
appear in Figures [31 and [U 

Hyperfinc splitting is the dominant broadening mech- 
anism for the europium spectral lines. The interaction 
between the nuclear spin and the atom's angular momen- 
tum vector causes energy level splitting in atoms with 
odd atomic numbers (europium is Z — 63). The effect 
is particularly pronounced in rare earth elements. The 
4129-A and 4205-A lines, for example, have FWHMs of 
1.5 A, due largely to hyperfine structure (but not isotope 




4129.0 4129.5 4130.0 

Wavelength (A) 



4130.5 



Fig. 2. — The 4129-A Eu II line in the solar spectrum. Individ- 
ual lines are annotated and are listed in Table [3] The hyperfinc 
components (see ti l 3 - 16 appear in the inset plot, which is aligned 
with the wavelength scale of the main plot. The re lative strengths 
of the 32 hyperfine components (Ivans et alJ 2006?) are plotted on 
a linear scale in the inset; the top half of the inset contains the 
components from the 151 Eu isotope while the 153 Eu isotope com- 
ponents appear on the bottom. The gray box indicates the portion 
of the spectrum used to calculate x 2 during the abundance fitting 
step (see § 13.21 1. The cross-hatched region indicates a portion of the 
spectrum used to fit a continuum. This plot represents a subset of 
the spectral region used in our analysis; the full region is 5 A wide 
and contains three additional continuum fitting regions. 




4204.5 4205.0 4205.5 

Wavelength (A) 



4206.0 



Fig. 3. — The 4205-A Eu II line in the solar spectrum. Individ- 
ual lines are annotated and are listed in Table [4] The hyperfinc 
components (see § 13.11 1 appear in the inset plot, which is aligned 
with the wavelength scale of the main plot. The re lative strengths 
of the 30 hyperfine components (Ivans et al. 2006) are plotted on 
a linear scale in the inset; the top half of the inset contains the 
components from the 151 Eu isotope while the 153 Eu isotope com- 
ponents appear on the bottom. The gray box indicates the portion 
of the spe ctru m used to calculate \ 2 during the abundance fitting 
step (see § 13.21 1. The cross-hatched region indicates a portion of the 
spectrum used to fit a continuum. This plot represents a subset of 
the spectral region used in our analysis; the full region is 8 A wide 
and contains six additional continuum fitting regions. 
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TABLE 3 

Line list for the region near Eu ii at 4129 A. 



A 

(A) 


hjlement 


Lower Level 
(eV) 


l°f 
solar fit 


VALD a 


r 6 

solar fit 


VALD b 


4129.147 


Pr II 


1.039 


-0.100 


-0.100 


-7.454 




4129.159 


Cr I 


3.013 


-1.948 


-1.948 


-6.964 


-7.362 


4129.159 


Ti ii 


1.893 


-2.300 


-1.730 


-6.900 


-7.908 


4129.166 


Ti I 


2.318 


-0.200 


-0.231 


-6.900 


-7.572 


4129.174 


Ce ii 


0.740 


-3.000 


-0.901 


-7.493 




i -i nn "loop 

4129.182 


Cr I 


2.914 


—0.100 




—8.300 




4129.220 


Fc I 


3.417 


—3.500 


—2.030 


—6.857 


—7.255 


4129.220 


Sm ii 


0.248 


— 1.123 


— 1.123 


—7.536 




4129.425 


Dy II 


0.538 


—0.522 


—0.522 


—7.554 




4129.426 


Nb I 


0.086 


—0.780 


—0.780 


—7.462 




4129.458 


Fc I 


3.396 


— 1.950 


— 1.970 


—6.863 


—7.206 


4129. 522 d 


Fe I 


3.140 


-3.497 




-6.873 




4129.643 


Ti I 


2.239 


-1.987 


-1.987 


-7.529 


-7.529 


4129.705 


Eu II 


0.000 


+0.260 


+0.173 


-7.174 




4129.817 


Co I 


3.812 


-1.808 


-1.808 


-6.099 


-7.782 


4129.837 


Nd ii 


2.024 


-0.543 


-0.543 


-6.237 




4129.959 d 


Fe I 


2.670 


-3.139 




-7.322 




4130.035 


Fe I 


1.557 


-3.900 


-4.345 


-7.885 


-7.826 


4130.036 


Fe I 


3.111 


-2.350 


-2.636 


-8.026 


-7.857 


4130.073 


Cr I 


2.914 


-1.971 


-1.971 


-6.929 


-7.349 


4130.122° 


V I 


1.218 


-1.000 


-3.142 


-7.060 


-7.800 


4130.233 


Mn I 


2.920 


-2.400 


-3.309 


-7.900 


-7.784 


4130.315 


V I 


2.269 


-0.300 


-0.607 


-7.187 


-7.585 


4130.364 


Gd ii 


0.731 


+0.177 


-0.090 


-6.608 




4130.452 


Cr 1 


2.914 


-1.099 


-2.751 


-6.805 


-7.348 



a All line data except Fq from Kurucz databases via VALD (unless other- 
wise noted). b Van der Waals parameters where available from Barklcm ct al. 
2000 via VALD (unless otherwise noted). c Identification from IMoore et aLl 
1966, log.qf and Tg from x 2 minimization. 11 Artificial iron lines included af- 
ter del Pcloso ct al. 2005.° Identification and parameters from Kurucz line lists 
hosted by the University of Hannover: 

http : //www.pmp .uni-hannover . de/ cgi-bin/ssi/test/kurucz/sekur .html. 

splitting — see insets of Figures [2] [4j . Since the relative 
strengths of the hyperfine components are constant with- 
out regard to temperature, pressure, or magnetic field 
(|AbtJll952l ). the components measured in laboratory set- 
tings can be applied to stellar spectra. 

Like other spectral fitting packages, SME has no built- 
in treatment of hyperfine structure. We therefore con- 
vert a single europium line into its constituent hyperfine 
components and include them as separate entries in the 
star's line list. The relative strengths and wa velength off- 
sets o f the hyperfine components come from llvans et all 
120061 which bases the values on an FTS laboratory a nal- 
ysis. Following the procedure of llvans et al.l l2006t we 
assume a solar system composition for the relative abun- 
dance of the two euro pium isotopes ( 151 Eu at 4 7.8% and 
153 Eu at 52.2%, from iRosman fc Tavlonll99l . We di- 
vide the Eu II log gf values listed in Tables 02 [H and 
[5] amongst the hyperfine components according to their 
relative strengths. All other attributes of the Eu n lines 
remain the same in the creation of the hyperfine lines. 
The relative strengths of the hyperfine components are 
displayed in the inset plots in the solar spectrum Figures 
El and 01 

3.2. Europium Abundances 

In order to measure the europium abundances in our 
selected stars, we begin with the co-added spectra de- 
scribed in § 12.21 We do a preliminary continuum fit us- 
ing the SME routines, which follow the VF05 procedure: 
deep features are filled in with a median value from neigh- 
boring spectral orders, then a sixth-order polynomial is 




6644.8 6645.0 6645.2 6645.4 6645.6 
Wavelength (A) 

Fig. 4. — The 6645-A Eu II line in the solar spectrum. Note that 
the ordinate axis is scaled differently than in Figurcs[2]and[3] Indi- 
vidual lines are an nota ted and are listed in Table [5] The hyperfine 
components (see §(3+J appear in the inset plot, which is aligned 
with the wavelength scale of the main plot. The re lative strengths 
of the 30 hyperfine components (Ivans ct al. 2006) are plotted on 
a linear scale in the inset; the top half of the inset contains the 
components from the 151 Eu isotope while the 153 Eu isotope com- 
ponents appear on the bottom. The gray box indicates the portion 
of the spec trum used to calculate \ 2 during the abundance fitting 
step (see C I3.2| |. The two cross-hatched regions indicate the por- 
tion of the spectrum used to fit a continuum. This plot represents 
a subset of the spectral region used in our analysis; the full region 
is 5 A wide and contains four additional continuum fitting regions. 
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TABLE 4 

Line list for the region near Eu ii at 4205 A. 



A 

(A) 


hjlcmcnt 


Lower Level 
(eV) 


l°f 
solar fit 


;gf 

VALD a 


r 6 

solar fit 


VALD b 


4204.695 


Y ii 


0.000 


-1.800 


-1.760 


-7.700 




4204.717 


Ce 11 


0.792 


-0.963 


-0.963 


-7.000 




4204.759 c 


CH 


0.519 


-1.140 


-1.158 


-8.900 




4204.771 c 


CH 


0.520 


-1.900 


-1.135 


-8.500 




4204.801 


Sm II 


0.378 


-1.771 


-1.771 


-6.738 




4204.831 c 


CH 


0.520 


-3.360 


-3.360 


-7.699 




4204.858 


Gd ii 


0.522 


-0.668 


-0.668 


-6.787 




4204.990 


Cr I 


4.616 


-1.457 


-1.457 


-6.658 


-7.852 


4205.000 


Fe I 


4.220 


-1.900 


-2.150 


-7.000 


-7.548 


4205.038 


V n 


1.686 


-1.850 


-1.875 


-6.800 


-7.913 


4205.042 


Eu ii 


0.000 


+0.250 


+0.120 


-6.800 




4205.084 


V ii 


2.036 


-1.100 


-1.300 


-6.900 


-7.956 


4205.098 


Fe I 


2.559 


-4.900 


-4.900 


-6.671 


-7.865 


4205.107 


Cr I 


4.532 


-1.160 


-1.160 


-6.582 


-7.776 


4205.163 


Ce II 


1.212 


-0.653 


-0.653 


-6.672 




4205.253 


Nd ii 


0.680 


-0.992 


-0.992 


-6.699 




4205.303 


Nb I 


0.049 


-0.850 


-0.850 


-6.677 




4205.381 


Mn ii 


1.809 


-3.300 


-3.376 


-6.800 


-8.001 


4205. 402 c 


CH 


0.488 


-2.300 


-3.960 


-8.000 




4205. 427 c 


CH 


1.019 


-2.300 


-1.130 


-8.000 




4205. 491 c 


CH 


1.019 


-3.463 


-3.463 


-8.000 




4205. 533 c 


CH 


1.019 


-1.800 


-1.149 


-8.000 




4205.538 


Fe I 


3.417 


-1.100 


-1.435 


-7.800 


-7.224 



a All line data except T§ from Kurucz databases via VALD (un- 
less otherwise noted ). b Van der Waals parameters where available from 
Barklcm et al. 2000 via VALD (unless otherwise noted). All molecu- 
lar line data (except Tq, from x 2 minimization) from Kurucz web site: 
http : //kurucz . harvard . edu/LINELISTS/LINESMOL/ . 



TABLE 5 

Line list for the region near Eu ii at 6645 A. 



A 

(A) 


Element 


Lower Level 
(eV) 


log 
solar fit 


gf 

VALD a 


r 6 

solar fit VALD b 


6644.320 c 


CN 


0.805 


-1.456 


-2.258 


-7.695 


6644.415 d 


La 1 


0.131 


-1.330 


-2.070 


-8.000 


6645.111 


Eu ii 


1.380 


+0.219 


+0.205 


-7.218 


6645.210 


Si i 


6.083 


-2.510 


-2.120 


-7.118 


6645.372 


Fe I 


4.386 


-2.759 


-3.536 


-6.780 -7.808 



a All line data except Te from Kurucz databases via VALD (un- 
less otherwise noted ). b Van der Waals parameters where available from 
Barklem ct al. 2000 via VALD (unless otherwise noted). c All molecu- 
lar line data (except Fq, from \' 2 minimization) from Kurucz web site: 
http : //kurucz .harvard. edu/LINELISTS/LINESMOL/. d Identification and pa- 
rameters from Kurucz line lists hosted by the University of Hannover: 
http : //www. pmp .uni-hannover .de/cgi-bin/ssi/test /kurucz/ sekur .html. 



fit to the region of interest. The built-in procedure cre- 
ates a flat, normalized continuum, though we find it nec- 
essary to fine-tune the continuum normalization in the 
course of our europium fitting. 

From the VF05 SPOCS catalog we take T oB , logg, 
[M/H], wsini, v mac , and v m i c (fixed at 0.85 km s _1 ) for 
each star we consider. In general, the global parameters 
from VF05 agree very well with the values adopted in the 
studies we compare to here. (Most of the literature val- 
ues fall within the 2-a errors quoted in VF05.) The VF05 
catalog is one of the largest and mos t reliable source s of 
stellar properties determined to date; lHavwoodll2006L for 
example, finds the VF05 T e s and [M/H] to be in good 
agreement with other reliable measurements. 

For most element abundance s we use a scaled solar 
system composition, shifting the lGrevesse fe Sauvallll99§ 



solar abundances by the star's [M/H]. The exceptions to 
this rule are sodium, silicon, titanium, iron, and nickel, 
which VF05 measured individually. Those a and iron- 
peak elements are therefore treated independently of the 
naive scaled-solar adjustment. It is possible that a more 
explicit treatment of iron-peak and a elements would 
improve our europium measurement accuracy. For the 
present study, however, we deem individual abundance 
analysis (apart from europium) unnecessary. We hold 
fixed the abundances of all elements other than europium 
in the subsequent analysis. 

We fit for the europium abundance by iterating three 
^-minimization routines that solve for the wavelength 
alignment, spectrum continuum, and europium abun- 
dance. A summary of each routine follows: 
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1. Wavelength. The pixel scale is pre-determined 
from the thorium lamp calibration taken each 
night. A first estimate of the rest frame wave- 
lengths comes from a cross-correlation of the full 
spectral segment with the solar spectrum, a built- 
in functionality of SME. We then use a 2- A region 
immediately surrounding the Eu II line of interest 
to perform a \ 2 minimization between the modeled 
stellar atmosphere and the spectral data, thus solv- 
ing for the wavelength scale alignment as precisely 
as possible in the Eu n region. 

2. Continuum. We fit a quadratic function across 
the points designated in the solar spectrum as 
continuum-fitting points (the cross-hatch regions 
in Figures [U [21 andH]). We adjust the quadratic 
continuum function vertically to require that 1-2% 
of spectral points in the full spectral segment are 
above unity, thereby ensuring that all spectra are 
scaled identically. 

3. Abundance. We perform a x 2 minimization ad- 
justing only the abundance of europium. We be- 
gin with the solar abundance value scaled by the 
star's metallicity, then search 1.0 dex of europium 
abundance space to find the best-fit value. In min- 
imizing the x 2 statistic, we calculate the residuals 
between the data and the fit in a limited region 
around the Eu II line (the gray regions in Figures 
IS andS]). 

The wavelength alignment, spectrum continuum, and eu- 
ropium abundance fitting routines are run in that order 
and iterated until a stable solution is reached. In most 
cases a stable solution requires only one or two iterations. 
The abundances for each line in each star as determined 
by this process are listed in Table [2] 

After running our automatic europium fitting algo- 
rithm on all 41 stars, we examine each line in each star 
by eye to confirm that the fit is successful. In a few of 
the metal-poor 4129-A fits, the blended lines that en- 
croach on Eu II were poorly enough fit with the VF05 
SPOCS values that the europium value was unconvinc- 
ing. In those cases, the 4129-A value is omitted from 
Table [2] and does not contribute to the average. See Fig- 
ure [5] for a comparison of a rejected 4129-A feature and 
a robust 4129-A fit. Similarly, in a few metal-poor stars 
the 6645-A line is too weak to be seen in the noise, and 
hence the output of the fitting routine serves only as an 
upper limit on the europium abundance. See Figure [5] 
for a comparison of a measureable 6645-A feature and a 
feature that provides an upper limit. In the cases where 
the 6645-A line is an upper limit only, it is listed as such 
in Tableland does not contribute to the average. 

If both 4129 A and 6645 A proved problematic we re- 
moved the star from our analysis entirely. The three stars 
for which this was the case (listed in £j |2.ip are omitted 
from Tables [T] and [2] Since all of the rejected stars were 
metal poor, we conclude that our fitting routine is most 
robust at solar metallicity, and becomes less reliable at 
[M/H] < — 1. Temperature may also play a role, as one 
of the rejected stars, HD 64090, has T cS = 7300 K; VF05 
determined SME to be reliable between 4800 and 6500 K. 

For each star we calculate a weighted average europium 
abundance value based on the three (or, in some cases, 



LO- 
CKS 
0.6 - 
0.4 ■ 

x 
3 

E0.2 

x> 

a> 

N 

B 




5 0.8 

z 

0.6 



Eu 



HD 103095 
[M/H] =-1.16 
|Eu/H] unreliable 
I i i i I i i i I i i i I 



0.4 1 

HD 22879 
0.2[ [M/H] =-0.76 
|Eu/H] = -0.55 

0.0 



1.00 



0.98 



0.96 



1.00 



0.98 



0.96 



0.94 




Eu ii 



HIP 57450 

[M/H] =-1.42 

[Eu/H] <-0.95 

' I ' ' ' I i ' ' I ' ' ' I 



HD 193901 
[M/H] =-1.19 
[Eu/H] =-0.98 



4129.4 4130.0 6644.8 

Wavelength (A) 



6645.4 



Fig. 5. — Bad fits (top) versus good fits (bottom) in compara- 
ble stars. In the upper left panel, the SPOCS stellar properties 
fit the lines adjacent to Eu II poorly enough that [Eu/H] is unre- 
liable. That line is removed from further analysis. In the upper 
right panel, the 6645-A line is buried in the noise, meaning the fit 
represents an upper limit to the europium abundance. The abun- 
dance upper limit is noted in Table [5| but does not contribute to 
the overall [Eu/H] measurement in the star. The lower two panels 
include fits to stars with similar characteristics to the stars in the 
upper panels, but where the 4129-A and 6645-A fits were more 
successful. 

two) spectral lines. Weighting the average is important 
because for stars with relatively few observations (e.g., 
HD 156365 in Figure [J), the weak line at 6645 A should 
be weighted significantly less than the more robust blue 
lines. Stars with a larger number of observations and 
higher S/N spectra (e.g., HD 144585 in Figure [TJ should 
have the 6645-A line weighted more strongly. In order to 
determine the relative weights of the three spectral lines, 
we tested the robustness of our fit by adding artificial 
noise to the spectra. We describe that process in § 13.31 

3.3. Error Analysis 

We begin our error analysis by adding to the data 
Gaussian-distributed random noise with a standard de- 
viation set by the photon noise at each pixel. We then 
fit the europium line again, using the same iterative \ 2 - 
minimization process described in § 13.21 and repeat the 
process 50 times. The standard deviation of the 50 Monte 
Carlo trials determines the relative weights of the lines in 
the average listed in Tabled with lower standard devia- 
tion lines (corresponding to a more robust fit) weighted 
more strongly. 

As expected, the results of the Monte Carlo trials show 
that the larger the photon noise in an observation, the 
less that line is to be weighted. We derive a linear relation 
between the Poisson uncertainty of an observation and 
its relative weight, based on the Monte Carlo results. 
Because the Monte Carlo trials are CPU intensive, we 
plan to use that relation to determine relative weights in 
the future. 

We estimate the uncertainty of our europium abun- 
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-1.0 -0.5 0.0 0.5 
[Eu/H], 4129 A 

Fig. 6. — A comparison of the Eu II values measured in the 
41 stars of this study, as listed in Table [3] Vesta is marked with a 
plus symbol. The solid line represents a 1:1 correlation; the dashed 
line is a best fit to the data. Points are omitted where the fit was 
poor in one of the Eu II li nes. Upper limits for the 6645-A line are 
marked with arrows. See S I4.1l for a discussion of the quality of the 
fit. 

dances to be the sum of the squares of the residuals for 
all three lines, where we assign the residuals the same 
relative weights we applied when calculating the aver- 
age. We also impose an error floor of 0.02 dex, added to 
the uncertainty in quadrature. 

The error floor comes from a comparison of abundance 
values from the individual lines (Figure [BJ discussed in 
more detail in ^ 14. 1 [) ; error bars of 0.03 dex on each mea- 
surement would make the reduced yx^ =1. A mini- 
mum error of 0.03 dex on each europium measurement 
translates into an uncertainty of 0.02 dex for the final 
averaged value. The 0.02 dex error floor is included in 
the uncertainty quoted in Table O 

As an initial test of the robustness of our fitting rou- 
tine, we measure the europium abundance in a spec- 
trum of Vesta, which serves as a solar proxy. The val- 
ues are listed in the last row of Table [2] The three 
Vesta europium abundances have a standard deviation of 
0.022 dex, indicating that systematic errors (^0.03 dex 
for Vesta, its offset from the solar value) may be a sig- 
nificant portion of the error budget. 

For the moment we absorb any systematic error with 
our random error estimates. The full sample of 1000 
stars, the analysis of which will follow this work, will 
allow a far more thorough investigation of the depen- 
dence of our results on various model parameters (T e g, 
logg 1 , etc.) than can be accomplished here. Therefore, 
we delay a substantive discussion of systematics until we 
have more europium abundance measurements in hand, 
though we touch on it again in § 14.11 and § 14.21 It is im- 
portant to note that since most of the CCPS stars are 
similar to the Sun, and we are treating each star identi- 
cally, our results will be internally consistent. 



3.4. Notes on Individual Lines 

3.4.1. Europium 4129 A 

The europium line at 4129 A is the result of a resonance 
transition of Eu II, and is a strong, relatively clean line. 
It provides the most reliable measurement of europium 
abundance in a star. Our fit to the Eu II line at 4129 A in 
the solar spectrum (described in § 13. 1| appears as Figure 
[2 with the 32 hyperfine com ponents (16 each from 151 Eu 
and 153 Eu; Hvans et ai]|2006f ) represented in the Figure^] 
inset. 

In the course of stellar fitting, the europium abundance 
is determined from the 4129 A line in all 41 stars ex- 
cept HD 103095. That star is quite metal poor ([M/H] = 
— 1.16) and cool (T e ff = 4950), factors that likely con- 
tribute to the poor fit in the lines adjacent to the Eu II 
line of interest (see Figure [5]). 

3.4.2. Europium 4205 A 

The Eu II line at 4205 A is the other fine-structure com- 
ponent of the resonance transition responsible for the line 
at 4129 A. Though it is as strong as the Eu n line at 
4129 A, contamination from embedded lines (see Figure 
[3]) has the potential make the 4205-A line less reliable. 
However, it is useful as a comparison line for the results 
from the 4129-A Eu n line. Our solar spectrum fit at 
4205 A appears in FigureEl with the 30 hyperfine compo- 
nents (15 each from 151 Eu and 153 Eu; Hvans et al.|[2006T) 
in the inset. Despite its blended nature, the fit appears 
sound in all 41 stars. 

3.4.3. Europium 6645 A 

The Eu II line at 6645 A is weaker than the lines in the 
blue, but it is also relatively unblended, making it worth- 
while to fit wherever possible. Our solar spectrum fit at 
6645 A appears in Figure SJ with the inset plot showing 
the 30 hyperfine compo nents (15 each from 151 Eu and 
153 Eu; llvans et al.ll2006l ). In HIP 57450, which is metal 
poor ([M/H] = —1.42) and has only one observation, the 
6645-A line was lost in the noise (see Figure [5|), and only 
the 4129-A and 4205-A lines contribute to the final eu- 
ropium abundance; the 6645-A line provides an upper 
limit only. 

4. RESULTS 
4.1. Comparison of Individual Lines 

We compare our results from the three Eu II lines 
in Figure [6] The Vesta abundances, included as a so- 
lar proxy, are consistent with the stellar abundances to 
within 0.03 dex. By assigning a measurement uncer- 
tainty of 0.03 dex to each line, the reduced \/x^ for each 
of the three plots is unity. This measurement uncertainty 
is the source of the error floor discussed in § 13.31 

Calculating the best-fit line for each comparison plot 
(the dashed lines in Figure [5]) lends insight into possible 
systematic trends in our analysis. In the blue, the two 
Eu II lines (4129 A and 4205 A) have no apparent linear 
systematic trend: the slope of the best-fit line is 1.03. 
The red Eu n line (6645 A), however, exhibits a minor 
systematic trend: the best-fit line in the red has a slope of 
1.10 relative to the blue, i.e., a 10% stretching of the blue 
abundance values about [Eu/H] = roughly reproduces 
the red abundance values. Since this systematic trend 
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Fig. 7. — A comparison of the final Eu II values from this work with literature measurements. The solid line represents a 1:1 correlation. 
The horizont al err or bars represent the uncertainties as quoted in the source literature. The vertical error bars come from o ur analysis as 
descr i bed in §EQ1 Two star s, H P 9562 and HP 228 79 (with [Eu/H] of +0.12 and -0.54, respectively), were included in the lBensbv et al.l 
2005, del Peloso et al. 2005, and Woolf et al. 1995 analyses, so those two stars are represented by three points along the abscissa at the 
same ordinate value. See Table [T] for the full list of the europium abundances plotted here, and see § 14,21 for a discussion of the quality of 
our fit based on this plot. 

is absorbed by the 0.03 dex error bars assigned to each 
point (necessary to make reduced \fx 2 unity even when 
comparing the non-systematic blue lines), we make no 
attempt to correct this systematic trend here. 

Between the relatively low measurement uncertainty 
needed to achieve reduced \fx 2 of unity and the minor 
systematic trend that only appears in the red, we con- 
clude that the europium abundance values derived from 
the Eu II lines at 4129 A, 4205 A, and 6645 A are consis- 
tent with one another. 

4.2. Comparison with Literature Europium 
Measurements 

In Figure [3 we compare our final europium abundance 
measurements to the literature values. The agreement is 
quite good. The literature values, plotted with the er- 
ror bars quoted in the original studies, appear as the 
abscissa. Our europium values, plotted with the error 
bars calculated in § 13.31 and listed in Table [H appear 
as the ordinate. The solid line represents a 1:1 corre- 
lation. Comparing the points to the 1:1 line, the reduced 



.. x 2 = 2.78. That va lue is dominated by HD 103095, the 
iSimmerer et al.l [20041 data point with the smallest error 

2 — 



bars, and omitting it makes the reduced yx 2 = 1.57. 

Adopting the global values (T c g, [M/H], logg) used 
by the comparison studies (instead of the VF05 values) 



drops the reduced y x 2 to 1.93, indicating that some 
of the scatter in Figure [7] is from the choice of stellar 
parameters. We emphasize that the VF05 stellar pa- 



rameters are reliable; we calculate reduced \J x 2 using 
the comparison studies' values in an attempt to separate 
how much of the disagreement in Figure [7] is from the 
europium abundance technique and how much is from 
the parameters adopted. 

Omitting the outlier mentioned at the beginning of this 
section, we examine the data in Figure[7]to search for sys- 
tematics in our results relative to the literature values. 
We consider the effects of a global offset in our europium 
values and a linear trend with europium abundance, find- 
ing that systematic offsets of — 0.1 dex and linear trends 
of — 40% are needed to make reduced a/x 2 " climb to 2. 
We conclude that our error bars have sufficiently charac- 
terized our uncertainties. 

Overall, we find reduced ^fx 2 to be dominated by 
the points at [Eu/H] < 0.5, consistent with our con- 
clusion in § 13.21 that our results are most reliable near 
solar metallicities , although the larger error bars on the 
Wool f~et al.l 1 19951 points make the correlation very for- 
giving near [Eu/H] = 0. Our spectra have very high 
S/N: -160 at 4200 A in a single observation, enhanced 
substantially by our co-adding procedure (§ 12. 2p . The 
automated SME synthesis treats all stars consistently, 
especially important for line blends in the crowded blue 
regions. For these reasons we believe that our europium 
abundance technique is accurate and robust and that our 
smaller error bars are warranted. 

We conclude that our abundance measurements are 
consistent with previous studies, not surprising since 
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most abundance techniques rely on the same Kurucz 
stellar atmosphere models. Because the majority of the 
points in Figures [7] fall near the 1:1 correlation, we also 
conclude that near [Eu/H] = our errors are ~ 0.03 dex, 
though at [Eu/H] < —0.5, the errors may be as high as 
0.1 dex. 

5. SUMMARY 

We have established that our method for measuring 
europium in solar-metallicity stars using SME is sound. 
The resolution and S /N of the Keck HIRES spectra are 
sufficiently high to fit the Eu n lines in question. The 
values obtained from the three europium lines are self- 
consistent, and our final averaged europium value for 
each of the 41 stars in this study are consistent with 
the literature values for those stars. 

By employing SME to calculate our synthetic spectra, 
we are self-consistently modeling all the lines in the re- 
gions of interest. Any blending from neighboring lines 
is treated consistently from star to star, adding robust- 
ness to our europium determination. Using SME has the 
added benefit of allowing us to adopt the stellar parame- 
ters from the SPOCS catalog. Our automated procedure 
ensures all stars are treated consistently. 

Having established a new method for measuring stellar 
europium abundances, we intend to apply our technique 
to 1000 F, G, and K stars from the Keck CCPS survey. 
Our analysis of europium in these stars will represent 
the largest and most consistent set of europium measure- 
ments in solar-metallicity stars to date, and will provide 
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